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Abstract 

we study the statistical properties of inhomogenous Poisson networks, we perform 
a detailed analysis of the statistical properties of Poisson networks and show that the 
topological properties of random cellular structures, can be derived from these models 
of random networks, we study both two and three dimensional networks with non uni- 
form density and show that with a class of symmetric distribution P(Ai , A2, ■ ■ • , A^v) 
Lewis and Aboav-Weaire laws are obeyed in these networks. 
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1 Introduction 



In the past three decades there has been a considerable interest in studying a class 
of non-equilibrium systems known collectively as Cellular Structures^ . The geometri- 
cal and dynamical properties of these systems are best displayed in the familiar pattern 
formed by a soap froth confined between two transparent plates. Other examples in- 
clude poly crystalline domains in metals, ceramics, magnetic domains, biological tissues 
and mono layers of fatty acids on surface of water. One can also mention other exam- 
ples from material science, like cracking in glazes, fracture and dewetting of polymer 
films above their glass transition temperature @-|H|- Despite the diversity in systems in 
which cellular structures are formed, numerous experiments have shown that the long 
time statistical behavior of these systems are characterized by certain universal, system 
independent laws. This means that topological and geometrical constrains influence the 
properties of these networks in a very essential way. For example, the fact that energy 
is associated with the length of the edges of cells, immediately leads to the conclusion 
that in the two dimensional structures all the vertices are 3-valent. Combination of 
this result with the Euler character formula shows that the average number of sides per 
cell, (i) is equal to 6: 

(t) = 6. (1) 

Similar considerations in 3-dimensions, where energy is associated with the area of faces 
of cells, proves that all vertices are 4-valent and that: 

(/) = tt^v (2) 
6 - (i) 

where (/) and (i) are the average number of faces per cell and the average number of 
edges per face respectively. Besides these properties which are a direct consequence 
of the Euler character formula, experiments have revealed a number of very general, 



properties among which the most important are : i) Von-Neumann P] , ii) Lewis [1C] 
and iii) Aboav-Weaire law [O]. Von-Neumann law refers to the rate of expansion of 
a single cell and is accounted for theoretically in a satisfactory way. It has also been 
generalized to curved surfaces where the curvature of the surfaces plays a role both in 
the dynamic of a single cell and in it's stability |12| |. The Aboav-Weaire law which is 
statistical in nature and refers to the correlation of adjacent cells has been shown to 
hold in random 3-valent graphs, viewed as planar Feynman graphs of </> 3 theory and 
solved by techniques of matrix models [13]. The authors of [14] have investigated the 
statistical properties of two dimensional random cellular systems in term of their shell 
structures. 

However there has been no explanation of the empirical Lewis law which states that 
the average area of cells has a linear relation with the number of cells, for large number 



of edges. However recently there have been some attempts [15] to derive both Lewis 
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Figure 1: 2D Poisson network, (a) Parallel columns in the x direction with width 
do. (b) These columns are divided into cells by edges through the points of independent 
uniform or non-uniform poisson distributions. The structure in c is used for obtaining 
the probability density of cover lengths, L c , for a given L.(d)Cover lengths of segment 
L in two distribution which partly overlap. 



and Aboav-weaire laws, from Poisson networks, i.e. networks based on a Poisson dis- 



tribution of horizontal segments between a fixed set of parallel lines. In ref.[15] it was 



shown that in such networks both the Aboav-Weaire and the Lewis laws were obeyed. 



The analysis of ref. [15] was however restricted to the uniform density case, i.e: where 
the distance between the horizontal lines were all equal. In this paper we extend this 
analysis to the non-uniform case and add one more ingredient of randomness to the 
network, where the distribution of vertical lines is also probabilistic. The paper con- 
sists of two parts, the first part (section 2), where we rederive the results of in a 
somewhat clearer way, and part two (section 3) where we extend these results to the 
non- uniform density case. 

2 Poisson Networks With Uniform Density 

The two dimensional Poisson networks are generated as follows (fig.l) |l5f l . One takes 
a family of parallel lines in the y-direction in the plane. The distance between lines does 
not affect the topological properties of the network, although it affects the geometry. 
This distance is taken as uniform and equal to do- Suppose there are N columns, 
C\, C2, ■ ■ ■ , Cat, in the network. The region between two successive columns is divided 
into cells. The division of each column C a , is based on a Poisson point distribution 
in the y-direction. An edge in the x-direction is taken through each Poisson point 
(P-point). In this way a 3-valent network is generated in the plane which, although 
different from the realistic cellular structures, is simple enough to be analyzed closely 
for the study of various statistical-topological properties of the network. 

The 3-dimensional tetravalent network are generated by a similar process (fig. 2). 
One takes an arbitrary 2-dimensional trivalent network(base network) on the xy plane 
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Figure 2: Three dimensional Poisson network, (a) an arbitrary planar network. (b)The 
columnar structure of three dimensional Poisson network. The cell in c has i = 4, 
F = 10, and height L. 

(fig. 2. a) and takes vertical planes (parallel to the z-axis) through each edge which 
divides the 3-dimensional space into prismatic columns. In each column one considers 
a Poisson distribution of uniform density in the z-direction and divides the columns into 
cells by planes perpendicular to the z-axis through each P-point (fig.2.b). The height 
of a cell, L, is the distance between adjacent P-points in the associated distribution 
(fig.2.c). The average cell height is unity. 

Consider one of the columns, say C a . The Poisson distribution of n points in a segment 
of length L is given by : 

P(n) = ^e' L . (3) 

Consider a particular set of columns C a , Cp, . . . , C 7 . The probability that in a given 
distance L, there are n a , np, . . . , n 7 segments respectively in columns C a , Cp, . . . , C 7 , 
is equal to 

$(L, n a , np,..., n 7 ) = P(n a )P(np) . . . P(n 7 ). (4) 

Now consider a reference column which we denote by C a . The probability that there is 
a cell with height between L and L + dL in this column is e~ L dL. I denote the columns 
which are neighbors of this reference column by C ai , C a2 , . . . , C ak . The number of 
these neighboring columns depend on the dimensionality and type of the network. The 
joint probability that there is a cell of height between L and L + dL in the column C a 
such that its neighbors C ai , C a2 , . . . , C ak respectively have n ai , n a2 , . . . , n ak points in 
the interval L is : 

, n ai , n a2 , • • • , flak 

)dL = P(n ai )P(n a2 ) . . . P(n ak )e- L dL. (5) 

Clearly this distribution function depends on the individual values of the variables 
n ai , n a2 , . . . , n ak , however in the sequel another probability distribution will be useful. 
That is ^(L, I)dL were I is defined as / = n ai + n a2 + . . . + n Qfe . The joint probability 
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fy(L, I)dL, is calculated as follows: 

* a (L, I)dL = ]T ^a(L, n ai ,n a2 ,..., n Qk )dL = ^-e~^ L dL : (6) 

where Yl' means the sum over all n Qi 's subject to the constraint that their sum be 
equal to I. It's important to note that the distribution ^f a (L,I) is obtained from a 
sum of distribution of the form ^ a (L , n ai , n a2 , . . . ,n a ) and is not equal to any of 
them, (compare with eq.(9a) of O] ) . 

2.1 Statistical Properties of Poisson Networks 
2.1.1 Two Dimensional Networks 

In two dimensions, every column has two neighbors, hence in eq.(||), k = 2, and the 
number of edges of a cell (see fig.l) in a reference column is equal to : 

i = / + 4. (7) 

Combining eqs.(||) and (|?|) gives the probability distribution g(L,i) of finding i-cells of 
height in the interval (L, L + dL): 

g(L,i)dL= \. ' e~ 3L dL. (8) 
(i-4)! 

Clearly this distribution function is normalized. 

/ Y.9(L,i)dL = l. (9) 

From eq.@ one can obtain the average height of i-cells: 

<L> 1= r L ^-dL, (10) 

Jo g(i) 

where g(i) is the total probability distribution of i — cells which given by: 

g®= I g(L,i)dL = -(-Y~\ (11) 



A simple computation now gives: 



o 



<L> 1= 1 -^. (12) 



This equation, then can be used to calculate the average area of i-cells. For the Poisson 
networks (fig.l), where the width of all cells are equal to do, we find: 

< A >,= d < L > t = , (is) 
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which as far as linear dependence on the number of edges is concerned, agrees with 
Lewis law. One can also define a generating function: 

oo 

G(9)=<e*>=X>(i)e* (14) 

i=4 

a simple calculation shows that: 

G ^ = J^> or lnG(q) = 4q-ln(3-2e«). (15) 

From which one obtains by successive differentiation various connected moments of the 
distribution: 

< i >= 6, (16) 

< i 2 > - < i > 2 = 6, (17) 

< i 3 > -3 < i 2 >< i > +2 < i > 3 = 30. (18) 

In the remaining part of this subsection we review briefly the basic steps of the analysis 
of ref.JO] for derivation of the Aboav - Weaire law. Consider an i-cell a of length L in a 
two dimensional Poisson networking. 3). This cell is in column and has two neighbors 
in this column, called a' and a". 

There are two adjacent neighbors 1 and 2, which respectively distribute n\ and n 2 
points inside the cell a. The number of sides of, o, is then equal to 

i = ni +n 2 + 4. (19) 

The total number of sides of the cells in column 1 adjacent to a is: 

J x = 4(ni + 1) + {Vi + 1) + (VI + 1) + V", (20) 

where the meaning of the numbers V\, V{ and V{' are specified in fig. 3. For obtaning the 
Aboav- Weaire law we use one of the three figures (fig.3a), (fig.36) or (fig.3c) and results 
are equivalent. Clearly inn which is the total number of sides of the cells adjacent to 
a is: 

inn = 12 + Ji + J 2 . (21) 

Where the definition of J 2 is similar to that of J\ and the number 12 comes from the 
average of sides of a' and a" . For the average value of V\, V{ and V" we use: 

< Vx + V{ >= <Lci ~ Lq> =< L C1 - L a >, (22) 

< V" >= < L * > =< L C1 - L a > + < L a > . (23) 

Where the segment L Cl in column Ci, which is the smallest segment containing L a , is 
called the covering length of L a [15]. With the probability distribution found in [15] 
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Figure 3: These figures are used for obtaining Aboav - Weaire law in two and three 
dimensions. 



for L c it is shown that the average < L c — L a >= 2, in any arbitrary column, it then 
follows that: 

< im, >= 24 + 4(ni + n 2 )+ < Vy + V{ > + < V 2 + Vi > + < V" > + < V 2 " > . (24) 
By substituting eqs.fflB), (E3), (E3) and <M) one finds that: 



<imi> = 16 + 4i + 2<L>j, 



16 + 4^ + ^^1 = 14 + — i. (25) 
3 3 



Combining eq.(25) with (12) one obtains 



< imi >= 16 + M + — — — = 14 + — i (26) 

3 3 



This is in accord with Aboav- Weaire law, i.e: < irrii >= c\ + c 2 i with c\ and C2 

L4 

3 



which are 14 and 4- respectively. 



2.1.2. Three-Dimensional Networks. 

Consider a 3-dimensional network and a column whose base is an i-cell. In eq. (|6j) , 
we should now equate k with i and, relate I with the number of faces as follows: 

I + 2 + i = F. (27) 

The probability distribution of F- cells (cells with F-faces) of height between L and 
L + dL, is now obtained as: 

9i(L,F)dL = ^f^- e -(H-i^dL. (28) 



The total distribution of F — cells whose base is an i — cell is given by: 

(» + l)^ + l 



oo 1 ? 

m (F) = I gi(L, F)dL = __(_-)*'-*- 2 . (29) 
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One can then obtain the average height of F- cells whose bases are i-cell. 



OO 



M^ = ^_l. (3 „) 
Jo 9i{F) + 

The average height of F-cells with any base is given by: 

oo 

< L > F = J2 <L >F,i 9(i), ( 31 ) 

i=4 

where g(i) is the probability distribution of i- cells in the base. As an approximation 
one can assume that the base network is hexagonal and obtain: 

F — 7 

< L > 6 ,f=< L > F = ——. (32) 

So the average volume (Lewis law in three dimensions) of F-cells will than be: 

F — 7 

<V > F =S——, (33) 

where S is the average area of cells in the base. Also, as second example, we assume 
that the network of base is Poisson network with uniform density, which has explained 
in previous section(eq,(|TT|)), so that eq.(31) will become: 



, oo 2 
< L >F = i E < L >F,i (o)^ 4 = 0-!55F - 1. (34) 

1=4 

Therefore by making use of eq.(^) we can obtain the overage volume of F — cells in 
this model: 

<v> J .= ^(i-3)(0.155F-l). (35) 
3 

From the formula for g%(F) various moments can be calculated: 

< F >i = 2i + 2, (36) 

< F 2 >i = 5i 2 + 9i + 4, (37) 

and 

< F > = 2 < i > +2 = 14, (38) 

< F 2 > = 5 < i 2 > +9 < i > +4 = 268. (39) 

3 Poisson Networks With Non-Uniform Density 

A little modification of the previous formulas is necessary to study the properties 
of a network where the density of points (the average distance between points to be 
denoted by A, in the following) in different columns are different. Instead of eqs.(|3|) 
and (|j) we have: 

Pa(n) = ( A~ r ^! e ~^' (40) 
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^ {L,n 1 ,n 2 ,...,n k ) = Pi{n 1 )P 2 (n 2 ) ■ ■ ■ Pk(n k ), 

k 

= IKfr^-A (4i) 



a=l* a Ua - 

respectively. By defining *f>(L,I)dL as before, we obtain: 

V (L,I) = $(L,ni,n 2 ,...,nk)^-e~*SdL, 

ni+ri2+...+n k =I 



Aq > -Lgp 



(42) 



where £o = ^ + ^ + XJ + -- - + XT;- Now we study another probability density. We 
consider a cell of length L in which randomly located in the C n th column (L is Poisson 
segment in C n th column). The cover length, L Cn+1 , of L is defined as the smallest 
length in C n+ ith column such as it is contain length L(i.e. L Cn+1 > L)(fig(lc)). Here, 
we wish to obtain the probability density of the cover length, L Cn+1 , whereas the main 
distribution in the columns of network is inhomogenous Poisson. So this new probability 
destiny must satisfy the following conditions: 

1) In the constant length L n+ \ of L, there are m points of C n +ith non-uniform Poisson 
distribution, i.e. 

P n+1 (m,L n+1 ) = (^±l)-J- e -^T. (43) 
A n+ i ml 

2) There are segments / and /' in C n+ ith non-uniform density 



P(n + i)(l) = T—e Wi, (44) 

A n +1 
1 



P(n + i)(l') = t e (45) 



and 



L Cn+1 = L n+1 + l + l', (46) 

L = L n+1 + t + t', (47) 

t < l<L Cn+1 -L+t. (48) 

3) 

V n+ i(L Cn+1 ,L)dL Cn+1 =l. (49) 



/; 



' L 

where V n +i(L Cn+1 , L) is the non-uniform probability density of cover length L Cn+1 . 
Therefore, the probability density of cover length in C ra +ith non-uniform Poisson dis- 
tribution with above conditions is: 
P n +i(L Cn+1 ,L,m) = 

J P n+1 (l)P n+1 (l')P n+1 (m,L n+1 )5(L Cn+1 -L n+1 -l-l')dldl', (50) 
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where Q is a normalization factor. One can find probability density of cover length for 
any m: 

CO -i j 

V n+1 (L Cn+1 ,L) = £ V n+1 (L Cn+1 ,L,m) = —^(L Cn+1 - L)e~^ (L -+^ L) . (51) 

m=0 A ™+ 1 

We see that, this distribution is independent of L n +i and t. The average value of 
(L Cn+1 — L) for fixed L in column C n+ i is: 

<(L Cn+1 -L)>=2X n+1 . (52) 

Note that, it is independent of L. In the case which we study the 3D poisson network, 
for a given segment L, there are two cover lengths, L Ck and L Cj in two independent 
poisson distribution, in which are neighbour with each other (figl.d). It is possible that 
these two cover lengths partly overlap. The probability density of the non-overlap 
length \L Ck — L Cj j is exp( Cfc Afc Cj ), for C^th column. The average of the non-overlap 
length in C^th column is: 

< \L Ck -L c .\ >= A fc . (53) 
3.1 Statistical Properties of Non-Uniform Density 

3.1.1 Two Dimensional Networks 

In a two dimensional network every column say C n has two neighbours, C ra _i and 
C n +i. The number of sides of i-cell in column C n , is i = 4 + 1, where / = n n -\ + ?^n+i 
and n n _i+n n +i are the number of vertices contributed by the cells in adjacent columns. 
Then for the distribution of i-cells of length between (L and L + dL) in column C n , we 
obtain: 

where £ n = + ^ + x^jTT- Hence, the distribution of i-cells with any arbitrary 
length in column C n is: 

/•OO 1 1 

9n(i)= g n (L, l )dL = -—(l-—y-\ (55) 

clearly this distribution function is normalized. From eq.(48) and (49) we also obtain 
the average height of i — cells in column C n : 

Jo g n (i) in 

So that for the Poisson networks with non-uniform density in one direction {y direc- 
tion) (figl) Lewis law (average area of % — cells) in column C n is: 

< A > n>l = . (57) 
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One can now evaluate various moments of the distribution g n {i). In all the following 
calculations, we use the following formulas for geometric sum: 

OO -I 

£*' = — x , (58) 



fc=0 



5>* = Tj^r, (59) 

fc=0 {L X) 

VW = ^±4. (60) 
to (i-^ 



A simple calculation shows that: 
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<i>n = J2Wn{i) = 3 + \ n i n = 4 + A n (- + - ), (61) 

i=4 A n-1 A„, + i 



DC 



< >n = £ i 2 g n (i) = 2(A„e„) 2 + 5A„£ n + 9, (62) 

where < i > n is the average number of i-cells in column C n . The variance turns out to 
be equal to: 

<i 2 > n - <i>l={\nZn){\ntn-l)- (63) 

In order to find the moments < i > and < i 2 > in the whole lattice and not in a 
particular column C n , we replace the average over the cells in a network, by an average 
over the densities in an ensemble of networks. Thus from eqs.(|57|), ( |6l|) and (??) we 
have : 

<A»i = d (i-3)< — >, (64) 

1 1 

<i» = 4+ < A„(- + - ) >, (65) 

< i 2 » = 2 < (A n £n) 2 > +5 < A n £ n > +9. (66) 

Where the averages on the right hand side are performed with a suitable distribution 

P(Ai,A 2 ,...,Atv). 

For simplicity in the following we assume that this distribution, is a symmetric distri- 
bution, i.e: 

P(Ai,A 2 ,...,A 7V ) = P(A 1 >,A 2 ,,...,A i vO, (67) 

where the primed indices are a permutation of the unprimed ones. Clearly this minor 
restriction allows our conclusions to be valid for a very large class of probability distri- 
butions. Finally, we consider Aboav-Weaire law for non-uniform Poisson distribution. 
By using of eqs.(|2^),([23|) and (|52[), we obtain: 

< Vt + V{ > n+l = < Lc - + \~ La >n+1 = 2, (68) 
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s V" ^ - _ Lcn + 1 >n+1 j- _ 1/(1 >n+1 o j_ 1 / r ^ teen 

< >n+l— 7 1 : — 1 + t < i^a >n+l • loyj 

An+l A n +i A n+ i 



Combining eqs.(p4]), ( |68| ) and (69), we have: 

< L a >( n +l),i ^ 1 

By using of eqs.(|57|) and (ffO[) , obtain: 



< im, >= 16 + 4i + Q . (n+1) '' + t < L a > {n -i),i . (70) 

An.J-1 An 



<im, >= i6 + 4t + (*-3)(- \ + - \ ). (71) 

A n +l?n+l A n _i^ n _i 



At last we have: 



< i m% >= 16 + 4i + (i - 3) < - — + - — > . (72) 

An+l?n+l A n __i£; n _i 



Where the last averaging are performed by a given distribution in eq.(67). If we use of 
uniform distribution we have A n+ i£ n+ i = A n _i£ n _i = 3, so that: 

14 

< irm >= 14 + —i, (73) 



we see that this result is the same in eq. (|26j) 

3.1.2 Three-Dimensional Networks 

Consider a 3-dimensional laminated Poisson network, based on an arbitrary 2- 
dimensional network, and a particular i — cell is the base. I label the column with 
this base by Cj o- The number of faces in a cell in the column above this i-cell is 

F = i + 2 + J, (74) 

where J = ri\ + n<i + . . . + rij, and is the number of additional lateral faces which 
results from the cells in the adjacent A:— column (fig. 3). Let fi(J) be the fraction of the 
F — cells in column C^o- In formula (41), one should now replace k by i, and I by J to 
obtain: 

£F-i-2/V l_\F-i-2 

where £o = ^ + X|" + xj + -- - + x~- Here Ao is the density of P-point distribution in 
the column C^o and A a (a = 1, 2, . . . , i) are the density of P-point distribution in the 
adjacent columns. ^o(L, F,i)dL is the probability of finding an F-cell in column C^o, 
whose height is between L and L + dL. Integrating over L, one obtains 

v& (F,i) = -i-(l--i-) M , (76) 

Ao?0 Ao£o 

this is the probability of finding an F-cell in column C^o- The average height of F — cells 
with i — cell base in column Co is: 

/■«> L\Eo(F, L, i)dL F-t-1 

< ^ >o,i,F= / T 7 p x = t , (77) 

JO *0{-r,l) SO 
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hence, the average height of F — cells with any base in column Co 

oo 

< L > ,f= < L >o,i,F go(i). (78) 



i=A 



Here, we can calculate eq . (]78|) for three different base networks. Assume that the 
network of base is: 
1) hexagonal 



F — 7 

< L > ,f= —z—, (79) 

SO 



2) uniform Poisson network 



<L>o,F=\Y,(^ A )(\r\ (so) 

3 - =4 Co 3 



3) non-uniform Poisson network 



oo 



< L > 0iF = J2( F J - )9o(i), (81) 

note that the non-uniform Poisson distribution in horizontal columns (2D Poisson net- 
work) and vertical columns are independent. Therefore we can easily calculate the final 
averaging in eq.(|74|). One can now find various moments of this distribution: 

<F> ifi = Ao& + « + l, (82) 
<F 2 > ii0 = (i + l) 2 + (2i + l)(A £o) + 2(A £o) 2 . (83) 

By the same averaging procedure as in the two dimensional case one can obtain: 

< L = (F-i-1) < — >, (84) 

Co 

< F >j = < A £o > +» + 1, (85) 

< F 2 >i = (i + l) 2 + (2i + 1) < A £o > +2 < (AoCo) 2 > • (86) 

The later averages are performed over the density distributions eq.(67). Note that for 
the distribution with uniform density in vertical columns one has: Ao£o = i + 1 and the 
eqs.flU), d|) and (86) reduce to: 

F — i — 1 

<L», >F = — , (87) 
i + 1 



= 2i + 2, 

<F 2 »i = (i + l)(5i + 4), (89) 



which are in agreement with equations (pQ} , (pq), (p8[). 

In this section, we derive the 3D Aboav - Weaire law in non-uniform Poisson network. 
Consider a 3D cell a with F faces, height L a and a base with i edges in reference 
column, Cq- We want find the total number of faces in cells adjacent to a with fixed F 
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and %. This quantity is {Fmp)^^. The number lateral walls in column k is m,. So that 
the total number of faces, which obtain from P-points within L a in first - neighbour 
columns to a is: 



J = E J k = E K™* + 2 )K + 1) + (V* + 1) + (Vl + 1)], (90) 
fe=i fc=i 

and, the total number of faces, which obtain from P-point within second- neighbour 
columns to a is: 

irrii— 3i 

J '= E Vl = (i mi -3i)<Vl>, (91) 
fc=i 

where, (imj — 3i) is the number of vertical walls of the adjacent cells that contact 
neighbours, and also there are another faces which obtain from P-points in pairs of 
adjacent column k, j that are adjacent to each other. So the total number of extra 
faces of such adjacenceis is 

J" = 2 ± [(n k + 1) + <\ L c k -Lc j \> i (g2) 
k=i Xk 

where the number 2 is due to two distributions are overlap in each vertical wall. Clearly 
one can obtain (FmF)o,i as 

< Fm F > ,;= 28 + J + J' + J" . (93) 

Hence the number 28 comes from the average of faces of two cells adjacent to a in the 
same reference column, Cq. By substituting eqs. ( |9Q|) , ([H|), ( |9~2"1) in (|93]), and using 
of©, (H), (|D, (|D and (0), we have 

< Fm F > 0i4 = m,[(l + Uf + (2 - ^±^) 4 - 2] + [4 - f]F + + 20, (94) 

and at last, we have 

«Fm f »i = mi{[i + i < — >]F + [2- (i + I) < — >]i - 2} 

+ {4 - 3i < \- >}F + 3»(t + 1) < — > +20. (95) 
so Co 



Hence the double average are performed by a symmetric distribution such as (70). It 
is remarkable that these averaging are performed with respect to those distributions in 
vertical columns(for example z direction), and are independent of distributions in base 
network. So that if we use of non-uniform poisson network as base network, we have 
two averaging, which are independent. 
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4 Conclusion 



My main conclusion is that Lewis and Aboav- Weaire laws apply to Poisson net- 
works. I have shown that by more randomizing these networks; (i.e. in the sense 
of random density distribution in different columns) one can still obtain the above 
mentioned laws. In fact the most features of these networks have in common with real 
cellular structures. But the only property of real systems in which all the angles around 
any vertex are equal is obviously absent in these networks. So we see that these model 
satisfy Lewis and Aboav- Weaire laws with the coefficients which are slightly different 
from the results obtained in [1,9]. Therefore despite this modification the main fea- 
tures of this model is still unaffected and one can derive the statistical and topological 
properties of random cellular structures from this modified model. 
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